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Abstract. By modifying the Chen and Wu ansatz, we have investigated 
some Friedmann models in which A varies as p. In order to test the con- 
sistency of the models with observations, we study the angular size - red- 
shift relation for 256 ultracompact radio sources selected by Jackson and 
Dodgson. The angular sizes of these sources were determined by using very 
long-baseline interferometry in order to avoid any evolutionary effects. The 
models fit the data very well and require an accelerating universe with a pos- 
itive cosmo logical constant. Open, flat and closed models are almost equally 
probable, though the open model provides a comparatively better fit to the 
data. The models are found to have intermediate density and imply the ex- 
istence of dark matter, though not as much as in the canonical Einstein-de 
Sitter model. 

PACS numbers: 04.20. Jb, 98.80.Es 



1. Introduction 

The cosmological constant A, which was originally invoked by Einstein to 
obtain a static solution of his field equations but was subsequently rejected 
by him after the realization that the universe is expanding, has fallen in and 
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out of fashion several times. However recent observations by Perlmutter et 
al (1999) and Riess et al (1998) strongly favour a significant and positive A. 
Their findings arise from the study of more than 50 type la supernovae with 
redshifts in the range 0.16 < z < 0.83 and suggest Friedmann models with 
negative pressure-matter such as a cosmological constant, domain walls or 
cosmic strings (Vilenkin 1985, Garnavich et al 1998). The main conclusion 
of these works is that the expansion of the universe is accelerating. 

Although a cosmic acceleration can also be accounted for by invoking 
inhomogeneity (though at the cost of the Cosmological Principle) (Pascual- 
Sanchez 1999, Dabrowski 1999), postulating a A-dominated model solves a 
lot of problems at once. The cosmological constant supplies the "missing 
matter" required to make Q tot = 1 (as suggested by the inflationary models, 
though on the basis of little observational evidence). It modifies CDM by 
putting more power on large scales, as is compatible with the CMB anisotropy 
limits. It also removes the inconsistency between the age of the universe and 
that of the globular clusters for larger values of the Hubble parameter Hq 
(the subscript zero denotes the value of the quantity at the present epoch). 

However, even with this significant value of A, we obviously face the 
problem that the upper limit of A from observations is 120 orders of magni- 
tude below the value for the vacuum energy density predicted by quantum 
field theory (Weinberg 1989; Carroll, Press and Turner 1992). (It is custom- 
ary to associate a positive cosmological constant A with a vacuum density 
p v = A/8nG.) 

Different types of solutions have been proposed to solve this problem and 
these can be classified into two categories. The first one, advocated by parti- 
cle physicists, implements some kind of adjustment mechanism, for instance, 
a counter term in the Lagrangian which goes away with the effective cosmo- 
logical constant. However, as remarked by Zee (1985), apart from the fine 
tuning of the parameters required by this method, there is no known sym- 
metry which would guarantee the effective A-term being zero. The second 
approach, advocated by general relativists, considers A as a dynamic vari- 
able. The link between the two categories is generally considered by studying 
different classes of scalar fields (Ratra and Peeble 1988). 

The second approach, which is essentially phenomenological in nature, 
has been extensively investigated in the past few years. It argues that, due 
to the coupling of the dynamic degree of freedom with the matter fields 
of the universe, A relaxes to its present small value through the expansion 
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of the universe and the creation of photons. From this point of view, the 
cosmological constant is small because the universe is old. 

Several ansatzes have been proposed in which the A term decays with 
time (Gasperini 1987, Freese et al 1987, Ozer and Taha 1987, Gasperini 1988, 
Peebles and Ratra 1988, Chen and Wu 1990, Abdussattar and Vishwakarma 
1996, Gariel and Denmat 1999). Of special interest is the ansatz A oc S~ 2 
(S being the scale factor of the Robertson Walker metric) by Chen and 
Wu, which has been considered/modified by several authors (Abdel- Rahman 
1992; Carvalho, Lima and Waga 1992; Waga 1993; Silveira and Waga 1994). 
Through dimensional considerations, in the spirit of quantum cosmology, 
Chen and Wu argued that one can always write the vacuum energy density 
as Mpj times a dimensionless quantity where Mp\ = (hc/G) 1 ^ 2 is the Planck 
mass. They therefore considered 
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where lp\ = (Gh/c 3 ) 1 ^ 2 is the Planck length. They noted that if one estimates 
So by cto, then n < 1 or n > 3 would lead to either too big or too small 
values of A compared to the observed limit. Consequently one needs to put 
n = 2, which gives the right value. 

Carvalho et al (1992) realized that (1) is not the only possible dynamic 
law for A given by the Chen and Wu ansatz. One may, for example, also 
write 
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where tp\ and tu = H^ 1 are Planck and Hubble times respectively. With the 
same argument as stated above, one finds A oc H 2 . 

Thus the Chen and Wu ansatz was generalized by Carvalho et al by 
considering A = aS' 2 +/3H 2 and later on by Waga (1993) by considering A = 
aS~ 2 + f3H 2 + 7, where a, (3 and 7 are adjustable dimensionless parameters. 

In this paper, we modify the Chen and Wu ansatz further. Since a non- 
zero A might be thought of as producing significant non-gravitational long 
range forces in the evolution of the universe and consequently influencing 
the large-scale dynamics, it is natural to hope that the sought dynamical 
law for A might also depend upon some cosmological parameter related to 
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the material content in the universe. Moreover, there is still room for any 
dimensionless parameter in the Chen and Wu ansatz. One parameter of 
interest is the density parameter Q, which is the density of the matter in 
units of the critical density, i.e., Vl = p/ p c = d>itGp/3H 2 . We thus consider 
the ansatz 

A = n n H 2 , (3) 

where n is a new dimensionless cosmological parameter to be determined from 
the observations. We here do not consider the S' _2 -dependence because, in 
view of the present estimate of A being of order of Hq, an independence 
seems more natural if A is to relax to its present estimate due to the expansion 
of the universe. 

Equation (3), which can alternatively be written as A = (8nnG/3)p or 
— (n/3)Q, puts A on the same footing as the energy density p and this 

corroborates the Machian view. (Here JIa = Pv/p c = A/3H 2 denotes the 

energy density of vacuum in units of the critical density.) 

Incidentally we note that if one estimates the present "radius" of the 

universe by So ~ cto ~ cHq 1 , then the present gravitational attraction 

4ncGp /3H roughly balances the cosmic repulsion AqcHq 1 provided A m 

In this paper we investigate cosmological models based on the dynamical 
law (3), in the framework of general relativity, and study their observational 
implications. The paper is organized as follows: The basic equations describ- 
ing the models are presented in Section 2. In Section 3, we study the angular 
size - redshift relation in the models, this providing one of the classical tests 
of the observational cosmology. For this purpose, we consider the data set 
of 256 ultra compact radio sources (with redshifts in 16 bins in the range 
0.5 - 3.8) from Jackson and Dodgson (1997). Our results are discussed and 
compared with previous works in Section 4. 

2. The models and the field equations 

In order to introduce the ansatz (3) and compare our conclusions with earlier 
works, we here consider the Friedmann-Lemaitre-Robertson- Walker models. 
For a perfect fluid, the models are fully specified by the scale factor S(t) and 
the curvature index k & {— 1,0,1} of the spatial hypersurfaces t = constant. 
Denoting the derivative with respect to the cosmic time t by a dot and using 
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units with c = 1, the models are governed by 
the Raychaudhuri equation: 



S AttG . . A 
-- = — (p + 3p) -- (4) 

and the Friedmann equation: 

S 2 k 8nG A 

^ + ^ = ^P+3- (5) 

We also assume that p = wp is the equation of state of the perfect fluid 
matter. Causality then requires —1 < w < 1. With the help of (3), equation 
(4) reduces to 

S 4ttG / 2n\ 

implying that inflationary solutions (i.e., solutions with negative deceleration 
parameter q = —SS/S 2 ) are possible for n > 3(1 + 3u>)/2. We also note that 
the models are open (k < 0), flat (k = 0) or closed (k > 0) according as 
fi^3/(n + 3), respectively, provided n > —3. For n < —3, the non-static 
model is always open. 

The Einstein field equations, via the Bianchi identities, imply 

(Ri - l -Rg^ = = -SrrG (t^ - ^g^j (7) 
leading to the energy conservation equation 

< i + 3 < 1 + w) 4 + 8sr ' (8) 

By using equation (3), this leads to 

p = cs -9(i+w)/(n+3)^ n ^_ 3) C = constant > 0. (9) 

For n = —3, the dynamic model yields (1 +w)p = with S — t and k = — 1. 
This leads to the following two cases, (i) Milne's model: p = p = A = 0, (ii) 
w — — 1 with the only constraint on p and A coming from equation (3) leaving 
their functional form otherwise free. However, w — — 1 just corresponds to a 
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second cosmological constant and will not be considered further. With the 
help of (9), the Friedmann equation (5) can be integrated to give 
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dy, (10) 



where we have chosen the origin of time at S — Sj 0) with S = (occurring 
at the first time) in the models without big bang and at — in those with 
big bang. Equation (10) implies that n > — 3 for k — or 1. In the case 
= 1, S is restricted by 



5 (2n-3-9«,)/(n+3) > Q/g^n + 3)GC7. 
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The integral in equation (10) can easily be evaluated for a general value 
of n if = 0. For 7^ 0, it can still be evaluated for particular values of n 
and we describe one such solution below. However numerical integration is 
possible for any particular value of n. 

For k — 0, equation (10) yields 



, 2nGC 

3(1 +W)il- * 



(2/9)(n+3)/(l+w) 



(12) 



(n + 3) 

Thus the scale factor S evolves as S oc t( n + 3 )/ 6 in the radiation-dominated 
(RD) era and S oc £ 2 ( Tt + 3 )/ 9 in the matter-dominated (MD) era, which is 
a simple generalisation of the canonical Einstein-de Sitter solution. The 
modified expression for the energy density, 
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18(1 + w) 2 7iG 
in the present case, has the same time-variation as 
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in the Einsten-de Sitter model. The expression for A in the present case 
becomes 



A 



An{n + 3) 2 
27(1 + w) 2 ' 



(14) 
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which matches with the natural dimensions of A. 



Solutions for k ^ : The solutions with n — 1 are as follows. 
k = 1, w = 1/3: 

t = VaS - S 2 + - sm-\2aS - 1) - (15) 
where a = (327rGC)/9. 

k = 1, iu = 0: 

i = a 4 |— cos a; sin 7 x — - cos a; sin 5 rr — — cos x sin 3 x 
1 6 24 

— — cosxsinx + — x\, (16) 
16 16 J ' K ' 

with 

asin 2 x = S 1/4 . (17) 

k = —1, W = 1/3: 



t = VaS + S 2 - - cosh~ 1 (2a5 + 1). (18) 



k = —1, W = 0: 



t = a 4 {cos 8 £ sin 7 x cos 6 x sin 5 x + — cos 4 x sin 3 x — — cos 2 rr sin x 

1 6 24 16 



with 



+ ^ ln(sec x + tanx)}, (19) 



,tan 2 a; = ,S 1/4 . (20) 



One can also approximate the solutions for a general n if there is no 
significant variation in the value of Q compared to that in H during the 
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expansion of the universe. (This assumption is trivially satisfied in models 
with n = 3(1 + 3w)/2 and obviously in flat models also.) By using equation 
(3), the Raychaudhuri equation (4) can be written as 

H+ ({3 + 9w - 2n}^ + l^j H 2 = 0. (21) 

If the variation in Q during a time interval is not significant compared to 
that in H, one can approximate the solution of equation (21) as 

t « mH-\ (22) 

where m = 6/ [6 — (2n — 3 — 9w)Q]. This gives an approximate age of the 
universe as to ~ itiHq 1 . Hence to > Hq 1 for n > 1.5. The corresponding 
approximate evolution of S turns out to be 

S oc t m . (23) 

Although we shall not use these approximate solutions given by equations 
(22) and (23) in further analysis, we shall see in section 4 that the estimates 
of t calculated from these equations are very close to those obtained from 
equation (10). 

In order to compare the models with observations and determine the 
relative contributions of the different cosmological parameters at the present 
epoch, it is necessary to recast equations (4) and (5) in the following forms: 

2 [g + ^ao] = fio, (24) 



1 + tt k0 = tt + tt A0 . (25) 

Here = k/S 2 H 2 is defined as the curvature parameter. Now the age of 
the universe follows from equation (10) as 



#ofo = j£[(i + £) fi o (2n - 3)/(n+3) - n 

where = S/Sq and fa = Si/S . 
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dfa (26) 
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3. Comparison of the models with observations 

We now want to study the consistency of our models with observations. We 
know that as we look back from our position at r = and t — t to some 
object at a radial coordinate distance r±, we are also looking back to some 
time ti < t and some expansion factor Si = S(ti) < So. Note, however, that 
neither r±,ti, nor Si are directly measurable; what are measured are physical 
properties like redshift, proper or apparent angular sizes, velocities, luminosi- 
ties etc. Since the various cosmic distance measures depend sensitively on 
the parameters of the models, the physical properties of the distant objects 
are also influenced by these parameters. In particular, the dependence of 
the angular size of a standard measuring rod upon redshift z depends 
upon these parameters. For this reason, the 9 - z relation was proposed as a 
potential test for cosmological models by Hoyle (1959). The original expecta- 
tion that the test would be able to distinguish between the geometries of the 
various cosmological models, was not fulfilled because of the intrinsic scatter 
and evolution of the sources. However, Kellermann (1993) argued that the 
evolutionary effects could be controlled by choosing a sample of ultra com- 
pact radio sources, with angular sizes of the order of a few milliarcseconds 
measured by very long-base line interferometry (VLBI). These sources, being 
short-lived and deeply embedded inside the galactic nuclei, are expected to 
be free from evolution on a cosmological time scale and thus comprise a set of 
standard objects, at least in a statistical sense. He used a sample of 79 such 
sources and showed that a credible 6 - z relation emerged. His results were, 
however, limited to showing that the corresponding 6 - z diagram favoured 
an Einstein-de Sitter model (with Q = 1), thus implying the existence of 
a large amount of cold dark matter. Jackson and Dodgson (1996) extended 
this work by including a A term and showed that Kellermann's data could 
also fit a low density, highly decelerating model with large negative A. 

Later on, more extensive exercise was carried out on a bigger sample of 256 
ultra compact sources by Jackson and Dodgson (1997). They then concluded 
that the canonical inflationary cold dark matter model (O — 1? — 0) 
was ruled out at the 98.5 percent confidence level and that a low-density 
Friedmann model with either sign of A was favoured. 

The original data set for the ultra compact sources used by Jackson and 
Dodgson was compiled by Gurvits (1994). The full sample included 337 
sources, out of which Jackson and Dodgson selected the sources with z in 
the range 0.5 to 3.8 for reasons discussed in their paper. These sources, 256 
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in number, were binned into 16 redshift bins, each bin containing 16 sources. 
Recently the compilation of Jackson and Dodgson has been used by Banerjee 
and Narlikar (1999) in the quasi-steady state cosmology. We also use this 
data set to study the 6 - z relation in the present models and check its 
consistency with observations. For this purpose, we now derive the © - z 
relation in the models. 

If the observer at r = and t = t receives the light from a source at 
a radial distance r\ with redshift z, the Hoyle's formula gives the apparent 
angular size 6 of the source as 

e = ^±4 (27, 

where d is the known (or assumed) proper size of the source and the coordi- 
nate radius r\ is given by 

v(n) = / 28 
/v ; Js /(i+z) SS V ; 

with 

T)(ri) = sin -1 r\, k — 1 
= ri, k = 

= sinh _1 ri, k = — 1. (29) 
Using (3), (5) and (9), equation (28) reduces to 



v(n) = ^ I % [(i + 1) ^ 9/(n+3) - n w v 

where ip = S /S. We also note that 



-1/2 



dip, (30) 



S = W\7^' k *° ^ 
n V "fco 



and equations (3) and (25) give 



n=£-(l + n k0 )-3. (32) 
"o 
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It is clear from equations (27) and (29) - (32) that, once we fix flko and 
Qo, the theoretical - z relation can be completely worked out for given 
d and H . We assume H = 65 km/s/Mpc and d — 10 pc (this does not 
have any consequences as there is still quite uncertainty in the value of H ) 
and calculate the theoretical Q(z) at the mean bin redshifts for a range 
of parameters Qq and Q^o- Using the observed values of 0, and the same 
standard errors a, of the ith redshift bin as used by Jackson and Dodgson, 
we compute x 2 according to 



This has 16 degrees of freedom as there are no constraints on the parameter 
space. We have performed an extensive investigation in order to find best-fit 
parameters. For this purpose, we have tried to minimise x 2 with respect to 
f2fc and f2 and find the best-fitting solution as fi^o = —0.22 and f2 = 0.49, 
which is an open model. However, flat and closed models also fit the data 
quite satisfactorily. This has been shown through some examples and the 
results have been plotted in Figures 1 and 2. 




(33) 
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QkO = -0.22 




0.2 0.4 0.6 0.8 

Figure 1. x 2 plotted against f2 for 4 typical models. The minimum 
value of x 2 decreases for lower density models. 
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Figure 2. Best-fitting logO/logz curves for 4 models, to be compared 
with Jackson and Dodgson' data points. 



4. Interpretation and conclusion 

Figure 1 illustrates four typical cases of flko'- for each Qko, we have plotted 
the values of x 2 against Q . The 95 percent line (with 14 degrees of freedom) 
intersects each curve and leaves sufficient portions of the curves below it, 
which shows a reasonably good fit. The minimum value of x 2 decreases for 
lower density universes. Moreover, in accordance with the standard model, 
we have lower density for the open model. However, we also note that, 
contrary to the claims of Jackson and Dodgson (1997), the closing of the 
universe with f^Ao does not abolish the need to postulate nonbaryonic matter, 
though the density is not as high as in Kellermann's model. Although the 
test is not conclusive in determining the curvature signature k, leaving all 
three signs almost equally probable, the fit of the data with the open model 
is excellent, with a probability exceeding 50 percent! The corresponding 
probability for the other optimum case (f^o — 0.3) is about 30 percent 
which also represents a good fit. The other two cases lie in between. The 
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predicted age limit in these optimum cases is about {l-lA)^ 1 . The formal 
best-fitting values are shown in Table 1. 

Table 1: Best-fitting values of the different parameters for four typical 
cases of flko- 





^0 


^AO 


n 


x 2 


Hoto 


-0.22 


0.49 


0.29 


1.79 


13.34 


1.05 





0.58 


0.42 


2.17 


13.61 


1.15 


0.1 


0.62 


0.48 


2.32 


14.00 


1.22 


0.3 


0.70 


0.60 


2.61 


15.96 


1.41 



It is, therefore, clear from Table 1 that the data favour an accelerating 
universe with a positive A which is interesting in view of the results from the 
supernovae data. However, it may be noted that a decelerating universe also 
cannot be ruled out completely. We mention, for example, a model Q = 0.42 
and Qao = —0.02 with d = 8.9 pc which has a fairly good fit (x 2 = 15.39) 
to the data. However, this model is not interesting on the grounds of age 
considerations. It may be noted that if we use the approximate solution given 
by equations (22) and (23) to calculate H Q t , the estimates never differ by 
more than 4.5 percent from those shown in Table 1. 

In Figure 2, we have plotted the theoretical - z curves for the four 
models discussed above and compared them with the Jackson and Dodgson' 
data points. The figure shows that the models fit the data very well. 
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